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PREFACE 


This  document  reports  the  work  performed  by  McDonnell  Douglas  Astronautics 
Company  for  the  Naval  Surface  Weapons  Center,  White  Oak  Laboratory,  Silver 
Spring,  Maryland,  for  the  development  of  an  aerodynamic  code  for  heat  transfer 
analysis  of  a  missile  at  high  supersonic  speed  and  nose  pressure  distribution 
of  a  missile  with  a  slightly  blunted  nose. 

For  heat  transfer  calculations,  the  code  has  been  modified  to  generate 
temperature  history,  equilibrium  temperature  and  the  rate  of  heat  transfer 
at  the  locations  specified  by  the  output  of  the  NSWC  computer  code  PING  and 
BING  for  stress  analysis  by  NASTRAN.  For  a  blunted  nose,  the  code  provides 
the  pressure  distribution  at  the  specified  locations  also. 

Theoretical  background  for  the  methods  used  in  the  development  is  described. 
For  convenience,  only  a  minimum  of  inputs  have  been  added.  Test  cases  for 
verification  purposes  are  presented  together  with  sample  inputs  to  reduce  the 
effort  of  preparation. 


A  complete  listing  is  included. 
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SUMMARY 


This  report  represents  the  work  performed  under  contract  number  N60921-80-C-0259 
for  the  Naval  Surface  Weapons  Center,  White  Oak  Laboratory,  Silver  Spring, 
Maryland,  under  Naval  Sea  Systems  Command  Sponsorship  for  the  development  of 
an  aerodynamic  code  for  analyzing  aerodynamic  heating  of  a  missile  at  high 
supersonic  speeds  and  the  pressure  distribution  on  a  slightly  blunted  nose. 

For  analyzing  the  aerodynamic  heating  of  a  missile,  a  number  of  subroutines 
have  been  added  to  the  existing  MDAC  aerodynamic  code.  For  missiles  with  a 
slightly  blunted  nose,  additional  calculations  have  been  incorporated  to 
provide  the  pressure  distribution  on  the  nose  as  required. 
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Section  1 


INTRODUCTION 


In  1979,  a  computer  code  for  aerodynamic  load  analysis  for  a  complete  missile 
(Reference  1)  was  completed  by  McDonnell  Douglas  Astronautics  Company  (MDAC) 
under  Contract  N60921-78-C-0203  for  the  Naval  Surface  Weapons  Center,  White  Oak 
Laboratory,  Silver  Spring,  Maryland.  It  analyzes  the  complete  missile  or  its 
components  (wing  or  body)  and  provides  pressure  distributions  at  given 
locations,  as  specified  by  the  outputs  of  NSWC  computer  code  PING  (Reference  2) 
or  BING  (Reference  3)  on  the  missile  surface. 

For  missile  flight  at  high  supersonic  speeds  (2  £  Mro  <  6)  the  problem  of 
aerodynamic  heating  becomes  important  and  needs  to  be  analyzed.  As  an  example, 
for  an  insulated  flat  plate  aligned  parallel  to  the  flow,  the  surface  temperature 
can  reach  anywhere  from  1000°R  to  3000°R  for  a  speed  range  of  Mot  =  2  to  6. 

Thus,  the  induced  thermal  stress  can  seriously  affect  the  structural  integrity 
of  the  missile  and  failure  of  a  missile  component  could  result.  Consequently, 
development  of  capability  for  analyzing  the  aerodynamic  heating  and  the  skin 
temperature  of  the  missile  is  needed  to  complement  the  pressure  load 
calculation.  Furthermore,  with  the  development  of  the  code  based  on 
the  linearized,  small -disturbance  equations  of  motion,  the  inability  of  the 
code  to  treat  bodies  with  a  slightly  blunted  nose  also  needs  be  remedied. 

To  accomplish  these  tasks,  the  code  was  substantially  modified  and  new 
subroutines  added.  The  principal  effort  consists  of  calculating  additional 
flow  field  quantities,  i.e.,  Mach  number  and  temperature,  at  the  specified 
locations  and  the  pressure  in  the  nose  region.  After  much  manipulation,  the 
desired  results  of  the  computation,  e.g.,the  temperature  history,  equilibrium 
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temperature,  heat  flux  and  nose  pressure  can  be  produced  as  punched  output  in 
the  NASTRAN  format  or  saved  on  tape  (temperature  history  of  element)  as 
specif i ed. 

To  minimize  the  development  effort  and  for  the  user's  convenience,  the  structure 
of  the  original  code  was  retained  to  a  large  extent.  Consequently,  only  a  few 
additional  input  parameters  are  required  for  analyzing  the  aerodynamic  heating 
or  the  pressure  distribution  on  the  blunt  nose.  For  high  supersonic  cases 
(Mm  £2.5),  the  code  is  restricted  in  its  application  to  a  complete  missile 
(body  and  wing)  as  was  developed  previously. 
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Section  2 


AERODYNAMIC  HEATING  OF  A  MISSILE  AT  SUPERSONIC  SPEEDS 

At  supersonic  speeds,  aerodynamic  heat  transfer  can  heat  up  the  skin  of  an 
airplane  or  missile  to  temperatures  at  which  the  allowable  strength  of  the 
skin  is  reduced.  Several  methods  have  been  developed  by  different  authors 
to  account  for  aerodynamic  heating  effects  at  high  speeds.  Among  these  are 
Eckert  (References  4,  5  and  6),  Wilson  (Reference  7),  Van  Driest  (References  8, 

9  and  10)  and  Spalding  and  Chi  (Reference  11).  Out  of  these,  Eckert's  reference 
temperature  method  was  selected  for  incorporation  into  the  aerodynamic  computer 
code  because  of  its  reliability  in  predicting  experimental  heat  transfer  results 
and  its  compatibility  with  analyses  used  in  the  original  MDAC/NSWC  Aerodynamic 
Code. 

Using  Eckert's  reference  temperature  method,  net  aerodynamic  heat  transfer  to 
a  missile  skin  is  calculated  from 

=  hc  (tr  -  Tw)  -  eo  Tw4  0) 

The  heat  transfer  coefficient,  hc>  is  calculated  from 

hc  =  (Cf*/2)(Pr)"2/3  p*  U6  cp  gc  (2) 

=  4.851  Cf*  p*  U^ 


The  skin  friction  coefficient,  Cf*9  is  calculated  as  follows 


(2a) 


where 


Cf  *  (R  <  105) 

£  6  6 

(3a) 

Cf  *  (R.  >  106) 

Tt  e& 

(3b) 

V +  icu  ‘  V)(109io  V5) 

J6  t  J6  6 

(3c) 

(105  <  Re  <  106) 

6 


*/»  *V2 


Cf^*  =  .664/(Re*/N£) 


(4) 


and 


Cf  *  =  .370/(Tog10  (Re*/Nt)) 


2.584 


(5) 


Reference  values,  indicated  by  an  asterisk  (*}  are  evaluated  at  a  reference 
temperature,  which  is  given  by 


T*  *  <TS  *  V'2  *  -22r  *T  MS  T6 


(6) 


values  for  and  Nfc  are 


V 


3  body  of  revolution 
1  planar  body  (wing) 


(7) 


V 


2  body  of  revolution 
1  planar  body  (wing) 


(8) 


i 
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The  governing  differential  equation  from  which  the  skin  element  temperatures 
and  heat  fluxes  are  calculated  is 


pwCpb  ST  =  hr  (TR  ‘  V  ‘  coTW4  (9) 

The  skin  equilibrium  temperature  for  each  skin  element  is  calculated  from 
Equation  (9)  when  the  rate  of  change  of  temperature  goes  to  zero 

hc  <tr  -  v  ■  “  \\  ('0) 

The  temperature  and  heat  flux  for  each  skin  element  at  time  j  are  calculated  from 
Tj  *  Tj-1  +  Chc<TR  -  Tj-l>  -  4Vvt» 

and 

«j  *  hc<TR  -  V  -  “V 

respectively. 

Computation  intervals  are  spaced  closely  together  during  early  times  and 
farther  apart  at  later  times  so  that  precision  will  be  retained  during  periods 
of  rapid  temperature  rise  but  excessive  computation  time  will  not  be  required 
after  the  period  of  initial  temperature  rise.  Temperatures  and  heat  fluxes 
are  printed  at  each  tenth  computation.  Computation  and  print  intervals  are 
shown  below 


Time 

AtJ 

Atprint 

(sec) 

(sec) 

(sec) 

0-2 

0.05 

0.5 

2-6 

0.10 

1.0 

6-10 

0.20 

2.0 

10-50 

0.50 

5.0 
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Section  3 

AERODYNAMIC  PRESSURE  OF  A  MISSILE  WITH  A  SLIGHTLY  BLUNTED  NOSE 


The  MDAC/NSWC  Aerodynamic  Code  as  developed  and  modified  in  the  past  few  years, 
is  based  on  the  linearized,  small  disturbance  equations  of  motion.  Consequently, 
its  application  is  limited  to  sharp  nose,  streamlined  bodies.  However,  it  is 
recognized  that  for  some  missile  configurations,  a  hemispherical  nose  is 
required  for  operational  reasons.  The  aerodynamic  characteristics  of  blunt 
nosed  configurations  can  be  considerably  different  from  those  of  the  sharp 
nosed  geometry.  In  particular,  the  pressure  distribution  on  the  nose  will 
greatly  influence  the  aerodynamic  coefficients.  An  appropriate  assessment  of 
the  effect  of  the  nose  is  needed. 


To  provide  the  capability  for  analyzing  the  nose  pressure  and  thereby  extending 
the  ability  of  the  aerodynamic  code  to  bodies  with  a  blunted  nose,  the  code  was 
modified  based  on  the  following  approach. 


In  the  subsonic/transonic  case,  the  pressure  distribution  over  a  spherical  nose 
is  calculated  from  the  Janzen-Rayleigh  relation  (Reference  17)  for  M^  <  0.75 

cp  =  1  -  Jsin%  +  -J-  (1  +  yrg-sin^  +  ^5^-  sin4e)  (11) 

and  the  empirical  relation  suggested  in  Reference  18  for  0.75  2  M  Z  2.0. 


where 


S(e) '  W cos  9  *  F<9-  MJ 


c,™u  =  c<'  «. 


and 


-2.27 


F(6,  MJ  =  {0.78  cos  0  -  0.95  exp[-2.235  (M^  -  1 )]}sin  6  (14) 

For  supersonic  flows,  the  modified  Newtonian  method  as  suggested  by  Lees 
(Reference  12)  in  a  survey  by  Isaacson  and  Jones  (Reference  13)  is  used. 

Ve)  ■  W  8  <>5> 

W(£-')/(K2)  <,6> 

.  P±1  „  2]£r  [ —  vi! - (17) 

P»  1  2YM/  -  (y-1) 

In  both  cases,  these  expressions  are  valid  only  in  the  neighborhood  of  the  nose 
region.  Taking  into  consideration  the  basic  assumptions  of  the  analysis  used 
in  the  code,  the  application  of  the  nose  pressure  option  must  be  limited  to 
slender  bodies  with  a  slightly  blunted  nose.  More  precisely,  the  ratio  of  the 
nose  radius  to  the  body  radius  must  not  greatly  exceed  0.1. 
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Section  4 


COMPUTER  PROGRAM 

In  1979,  under  the  contract  number  N60921-78-C-0203  from  NSWC-WOl,  the  MDAC/NSWC 
Aerodynamic  Code  (Reference  1)  was  developed  for  calculating  the  aerodynamic 
pressure  distribution  over  the  missile  traveling  at  high  supersonic  speeds 
and  generalized  attitudes.  To  extend  further  the  capability  of  the  code,  it 
has  been  modified  to  include  the  aerodynamic  heating  at  supersonic  speeds. 
Additionally,  analysis  for  the  pressure  distribution  over  a  slightly  blunted 
nose  has  also  been  incorporated. 

As  before,  the  code  accepts  as  input  the  output  of  NSWC  computer  codes  PING  and 
BING  which  prescribe  the  finite  element  on  the  missile.  The  code  generates 
as  output  the  pressure  distribution  and  heat  transfer  results  for  subsequent 
stress  analysis  by  NASTRAN.  At  the  user's  option,  both  the  pressure  load  and  the 
equilibrium  temperature  may  be  produced  as  punched-card  output  in  the  NASTRAN 
format.  For  convenience,  results  af  the  aerodynamic  heating  analysis  are 
saved  on  tape  (TAPE  12)  for  future  studies. 

4.1  CODE  CAPABILITIES  AND  RESTRICTIONS 

Using  appropriate  input  values,  the  code  can  be  utilized  to  analyze  a  missile 
flying  under  a  combination  of  the  following  conditions: 

(a)  Circular  or  near  circular  body 

(b)  Single  or  multiple  wings 

(c)  Subsonic  or  supersonic  speed  (Mm  ?  6.0)* 

(d)  Small  angle  of  attack,  <10° 

(e)  Small  angle  of  yaw,  <_10° 

*  For  higher  subsonic  case  (M^  >  2.5),  code  applies  to  complete  missile  only. 


(f)  Sharp  or  slightly  blunted  nose,  r^/ rb  <  0.1 

(g)  Aerodynamic  heatinq  at  supersonic  speed 

Although  the  small  disturbance  basis  of  the  code  generally  excludes  operation 
in  the  transonic  regime,  e.g.,  0.8  <  Ma<  1.2,  the  blunt  nose  relations 
include  an  expression  for  pressure  coefficient  valid  in  this  regime  (Equation 
1?). 

As  a  result  of  the  basic  assumptions  used,  some  limitations  on  the  use  of  the 
code  need  be  noted.  For  a  missile  with  a  blunted  nose,  aerodynamic  heating  of 
the  nose  region  represents  an  extremely  complex  phenomenon.  The  present 
approach  selected  for  the  heat  transfer  anaysis  is  valid  only  for  analyzing 
the  sharp  nose  case  and  should  not  be  used  for  the  blunt  nose  case. 

The  inputs  specifying  the  missile  dimensions  and  the  free  stream  pressure  must 
be  consistent  with  each  other.  As  an  example,  when  the  missile  dimensions  are 
given  in  inches,  free  stream  pressure  (PINF)  must  be  in  pounds  per  sq.  in. 
However,  in  the  case  of  heat  transfer  analysis,  the  missile  dimension  must  be 
in  inches  and  the  free  stream  pressure  in  psi. 

Based  on  the  operating  experience  with  the  CDC  CYBER  174  computer  with  the  NOS 
system,  the  code  requires  a  field  length  of  166,000  for  execution.  The 
computing  time  varies  between  one  to  three  minutes  for  different  cases. 

4.2  INPUT  DESCRIPTION 

Several  input  parameters  have  been  added  to  the  input  list  (NAMELIST/UNIFID) 
for  specifying  either  the  heat  transfer  analysis  or  the  pressure  distribution 
analysis  for  a  slightly  blunted  nose.  The  original  input  parameters  retain 
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their  function  as  explained  in  great  detail  in  Reference  1  (pp.  9-16)  and  are 
repeated  here  for  convenience. 

XMACH  Mach  number  of  flow.  Higher  supersonic  routine 

employed  when  greater  than  2.5,  default  value  «  2.0. 

PINF  Free  stream  pressure,  default  value  =  14.7  psi. 

OAOEG  Incremental  angle  of  attack  of  missile  in  degrees, 

default  value  =  0. 

SIO  Case  identification  number  for  each  missile  component 

(up  to  20  allowed);  default  values  have  been  set  at 
1.  Also  used  for  selecting  the  sign  convention  for 
pressure  output.  Code  uses  the  accepted  rule  for 
pressure,  i.e.,  positive  acting  toward  the  surface  and 
negative  acting  away  from  the  surface.  For 
opposite-sign  convention,  set  SID(L)  as  negative  for 
pressure  output  on  the  desired  component  L  (in  the 
order  of  body,  wing  surfaces,  upper  or  lower). 

YAW  Angle  of  yaw  in  degrees  with  respect  to  the  flow. 

BROLL  Missile  roll  about  body  axis  (degree). 

The  following  input  parameters  are  required  for  the  case  where  missile  body  is 
included: 

NBPFL  Number  of  body  profile  stations  for  specifying  body 

geometry  inputs  r  =  r(x);  up  to  51  stations  allowed. 

X8  <  51  Axial  coordinate  of  body  stations  in  the  direction  of 

flow  for  specifying  the  body  profile. 

RB  _<  51  Radial  coordinate  of  body  profile  at  corresponding 

body  stations. 

ZD  <  51  Body  camber  at  corresponding  body  stations. 
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ARB  Angle  of  attack  of  missile  body  in  degrees,  default 

value  =  0. 

In  accordance  with  the  overall  organization  of  the  code,  the  input  parameters 
and  data  have  been  organized  into  two  main  groups  in  accordance  with  their 
assigned  functions.  The  first  group  contains  all  the  parameters  needed  for 
the  specification  of  case  selection,  missile  configuration,  and  flow 
conditions.  They  form  the  input  parameter  of  the  namelist  UNIFID.  The  second 
group  specifies  the  finite-element  and  grid-point  system  as  generated  by  the 
NSWC  computer  codes  BING  and/or  PING.  It  includes  the  specification  for  the 
coordinate  system  used  for  the  finite  elements. 

For  the  convenience  of  the  user,  card  format  for  inputs  are  described  in  their 
sequential  order  as  required  by  the  code. 

A.  Aerodynamic  Control  Parameters. 

Card  1  Title  and  case  description,  format  (20A4),  alphanumeric. 

Card  2  SUNIFID,  begins  at  Column  2  with  the  symbol  S. 

Card  3  And  additional  cards  if  necessary,  contain  the  following 
as  input: 

ICASE  =  1  Wing  only. 

=  2  Body  only. 

=  3  Complete  missile. 

N TRANS  Number  of  coordinate  transformations  required  to  convert 

the  coordinates  for  the  finite  element  to  the  coordin¬ 
ates  used  in  the  aerodynamic  code.  Two  transformations 
are  allowed  for  each  missile  component,  in  the  following 
order;  body,  wings,  and  nacelles. 
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I PUNCH  -  0 


No  punched  card  output,  default  value. 

«  1  Punched  card  output  requested. 

IRW  ■  2  Computation  results  are  not  saved  on  tape,  default  value. 

>  1  Aerodynamic  computation  results  are  saved  on  tape  1?  for 

restart,  first  run  only. 

*  2  Restart  run,  bypass  all  aerodynamic  computations  and 

begins  computation  at  the  start  of  pressure 
interpolation. 

POLAR  Number  of  incremental  angles  of  attack,  default  value  = 

0. 


The  following  inputs  are  required  for  the  case  where  the  missile  wing  is 
included: 


I SOLID  »  0  ,  Wing  of  built-up  construction,  pressure  load  on  both 

upper  and  lower  wing  surface  generated,  default  value. 


=  1 

I  FORM  =  0 


=  1 


*  2 

NWPl 


Wing  of  solid  section,  sum  of  the  pressure  loads  on  wing 
surface  generated. 

Wing  of  built-up  construction,  both  upper  and  lower 
surface  height  (ZWI)  are  given  at  identical  locations 
(XWI,  YWI) . 

Wing  of  built-up  construction,  upper  and  lower  surface 
height  (ZWI)  are  not  given  at  identical  locations. 

Flat  wing  surface. 

Number  of  coordinate  points  specify  the  wing  surface 
contour  (£30).  Code  allows  a  maximum  of  20  wing 
surfaces  (upper  and  lower).  Each  wing  surface  must  be 


NWING 

XWI(i,  j), 
YWI ( i ,  j), 
ZWI(i,  j) 
DIHEO 
PIVOT 

ARW 


defined  in  the  form  of  equation  of  a  surface,  ZWI  = 
f(XWI,  YWI).  A  maximum  of  30  points  are  permitted  with 
the  provision  that  the  first  four  points  are  restricted 
to  the  specification  of  the  corners  of  the  wing  only. 
The  coordinates  for  the  corners  must  be  given  in  the 
following  order.  Starting  with  the  innermost  point  of 
the  leading  edge,  the  remaining  corners  are  to  be  given 
in  a  counter  clockwise  order  viewing  from  the  above. 
Number  of  wings  or  tails.  A  maximum  of  10  is  allowed, 
default  value  =  0. 

Wing-surface  coordinates  where  indices  i  and  j 
designate  the  surface  coordinate  point  and  wing 
surface,  respectively. 

Dihedral  angle,  degrees. 

Indicates  dihedral,  0  for  no  dihedral,  and  >0  for 
dihedral  present. 

Angle  of  attack  of  wing  in  degrees,  with  respect  to  the 
body  axis  if  complete  missile  is  considered,  default 
value  =  0. 


The  following  input  parameters  are  required  for  the  complete-missile  case  (ICASE 
=  3)  except  for  high-supersonic-speed  (M  >  2.5)  case: 

NWPANL  Number  of  wings  directly  connected  to  the  body. 

XLG  Axial  distance  at  the  Intersection  of  wing  leading  edge 

and  body  for  each  wing. 


13 


XT’S  Axial  distance  at  the  intersection  of  wing  trailing  edge 

and  body  for  each  wing. 

Card  4  SEND  Begin  at  column  2  the  symbol  S  for  ending  input. 

8.  Finite-Element  and  Grid  Inputs.  The  second  group  of  input  cards  are 
arranged  immediate  following  Card  4,  as  follows: 

Cards  5  and  6  Coordinate-system  specifications  in  the  format  of 

NASTRAN  bulk-data  deck  [page  2.4-49  -  2.4-54,  NASTRAN 
User's  Manual,  NAS  SP-221(01),  1972].  Three  position 
vectors  -A,  B,  and  C  -  are  used  to  define  the  coordinate 
system.  The  first  defines  the  origin,  the  second 
defines  the  Z  axis,  and  the  third  defines  a  point  in  the 
XZ  plane. 

For  Card  5  the  following  format  is  used: 

Col  1-8  Coordinate  system  CORDJRbb  for 

rectangular,  CORDJCbb  for  cylindrical, 
CORDJSbb  for  spherical,  where  J  identifies 
the  coordinate  system  numbers. 

Col  9-16  Coordinate  identification  number  J,  integer. 
Col  17-24  Reference  coordinate  system,  integer, 
optional 

Col  25-32,  33-40,  41-48,  Components  of  vector  A 

(3F8.2). 

Col  49-56,  57-64,  65-72,  Components  of  vector  B 

(3F8.2). 
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For  Card  6:  Col  1-8  blank. 

Col  9-16,  17-24,  25-32,  Components  of  vector  C 

(3F8.2) . 

Cards  5  and  6  are  to  be  repeated  as  many  times  as  there  are  coordinate  systems 
(NTRANS)  required  for  specifying  the  finite  element.  They  are  to  be  arranged 
in  the  order  of  body  and  wings. 

The  finite  element  and  grid  specifications  are  arranged  immediately  following 
the  coordinate-system  input  in  the  same  order,  i.e.,  body,  wings. 

Card  7:  Specifies  the  number  of  grid  cards  (NGRIDP)  and 

element  cards  ( NE ) ,  Format  (2110). 

The  grid  and  the  element  cards  as  generated  by  the  NSWC  computer-code  BING  are 
sorted  into  two  separate  groups.  They  are  then  placed  immediately  after  card  7 
with  the  grid  cards  in  front. 

Some  remarks  regarding  the  input  parparation  need  to  be  made  here  to 
facilitate  the  use  of  the  aero  code. 

For  the  body-only  case,  inputs  of  the  body  profiles,  expressed  as  r  *  r(x), 
for  aero-panel  generation  must  be  in  the  aero-coordinate  system.  It  is 
required  that  x  be  the  axial  distance  from  the  tip  and  r  the  radius  of  the 
body  cross  section  (Figure  4-1).  For  the  finite  element  and  grid  point,  their 


^  •  Axit  in  »h»  Flow  Direction 

-  Axil  to  tha  Right 
Z%  ■  Axit  Politiva  Upward 

Figure  4-1.  Body  Profit  Input 

relationship  to  the  basic  coordinate  system  and  the  basic  coordinate  system  to 
aero-coordinate  system  are  to  be  specified  in  the  coordinate  input  card 
immediately  following  t:.e  namelist  input.  The  grid  card  and  the  element  card 
are  to  be  separated  and  placed  after  the  coordinate  system  specification,  with 
a  card  specifying  their  numbers  placed  behind  the  coordinate  card. 

The  input  for  the  wing-only  case  requires  the  specification  of  both  the  wing 
planforin  and  its  surface  contour.  The  input  for  the  description  of  the  wing 
planform  requires  the  coordinates  of  corner  points  and  break  points.  For  each 
individual  wing,  inputs  for  those  of  quadrilateral  shape  are  their  four  corner 
coordinates.  The  codes  have  allocated  the  first  four  values  of  the 
coordinates  XWI,  YWI  and  ZWI  for  each  wing  or  wing  region,  arranged  in  a 
counterclockwise  order  starting  from  the  leading-edge  corner  of  break  point. 
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XWII4.2) 

YWM4.2) 

ZWII4.2) 


XWI0.1) 

YWId.1l 

zwid.i) 


XWII2.1) 

YWII2.il 

ZWK2.1I 


XWK4.1) 

YWK4.1) 

ZWH4.1) 


XWK3.1) 

YWK3.1I 

ZWK3.1I 


XWII1.2) 

YWIM.2) 

2WK1.2I 


XWH2.2)  XWI  (3,2) 

YWM2.2)  YWI(3,2) 

2WI(2,2)  ZWK3.2) 


YWI(1.2> 

VWK1.3) 

ZWtfl.3) 

k 

XWK2.3)  XWM3.3).  XWK4.3) 
YWII24)  YWII3.3),  YWII4.3) 
ZWK2.3)  ZWK3.3),  ZWK4.3) 


NWING  -  3,  NTRANSU)  *  1.  1,  1, 

If  Built-Up  Wing  it  U«d,  NWPKtl,  NWPK2I.  and  NWPK3)  Mutt  be  Specified 


Figure  4-2.  Subdivifion  of  Complex  Wing 


In  the  case  of  triangular  wings,  there  are  only  three  corner  points.  To 
satisfy  the  same  requirement,  the  apex  or  wing  tip  is  to  be  considered  as  two 
coincident  points  of  the  same  coordinate.  The  inputs  are  to  be  arranged  like 
those  of  the  quadrilateral  shape,  with  leading-edge  body-corner  point  as  the 
starting  point.  Wings  with  complex  shapes  (see  Figure  4-?)  must  be  subdivided 
into  appropriate  quadrilateral  or  triangular  regions  and  the  inputs  prepared 
accordingly.  The  only  requirement  is  that  the  control  chords,  apart  from 
either  the  leading  edge  or  the  trailing  edge,  must  be  in  the  streamwise 
direction. 


All  inputs  for  the  wing  surface  are  to  be  given  in  the  form  of  the  equation  of 
surface  Z  =  f(x,  y)  in  terms  of  XWI,  YWI  and  ZWI.  Starting  with  the  fifth 
value,  XWI  (5,1),  YWI  (5,1),  ZWI  (5,1),  their  total  number  may  not  exceed  26 
points.  For  wings  with  flat  surfaces,  only  the  first  four  points  are 
necessary. 
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Inputs  defining  wing  thickness,  camber  and  twist  have  been  eliminated  to 
reduce  the  input  preparation  effort.  The  code  generates  the  required  values 
of  camber,  thickness  slope,  and  twist  angle  for  the  wing  using  the  existing 
surface-fit  scheme  from  the  input. 

Similar  to  the  body-only  case,  the  finite-element  and  grid-point 
coordinate-system  cards  are  placed  immediately  behind  the  namelist  input. 

They  are  followed  by  a  card  specifying  the  number  of  grid  cards  and  element 
cards.  The  grid-point  and  element  card  as  generated  by  PING  occupies  the  last 
position  in  the  input. 

The  coordinate  system  used  for  the  winq  can  be  any  one  of  three  systems,  i.e., 
rectilinear  (x,y,z),  cylindrical  (r,e,z)  and  spherical  (r,e,|i).  In  terms  of 
the  input  parameters,  they  are  represented  by  XWI,  YWI  and  ZWI,  respectively, 
for  all  three  systems. 

For  the  complete-missile  case,  the  input  can  be  considered  as  a  combination  of 
these  two  cases,  i.e.,  the  body-only  and  the  wing-only  cases.  Only  the 
parameters  ICASE  and  NTRANS  need  be  corrected  to  reflect  the  changes. 

Both  coordinate-system  cards  and  finite-element  cards  are  to  be  grouped 
together  in  the  order  of  body,  wings,  and  nacel’es. 

As  a  result  of  the  approach  used  for  the  analysis,  the  computations  for  the 
hiqh-supersonic-f low  case  are  considerably  simpler.  Accordingly,  the  code 
treats  the  complete  missile  as  a  rule. 
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As  previously  noted,  for  aerodynamic  heating  analysis,  inputs  for  missile 
dimensions  such  as  XB,  RB,  ZD,  XWI,  YWI,  ZWI,  XLG,  XTG  and  the  finite  element 
grid  inputs  must  be  given  in  inches.  The  free  stream  pressure  PINF  should  be 
in  psi. 

The  additional  inputs,  required  for  either  the  heating  transfer  analysis  or 
the  blunted  nose  case  are  described  below. 

4.2.1  Aerodynamic  Heating  Analysis 


IHTRANF 

=  0 

No  heat  transfer  analysis,  default  value 

=  1 

Heat  transfer  analysis  requested 

TI 

=  518. 0°R 

Free  stream  temperature,  default  value 

TEMPI 

=  540. 0°R 

Initial  temperature  of  missile,  default  value 

BMATL 

Body  skin  material  flag 

=  1.0 

Aluminum  (At), default  value 

=  2.0 

Titanium 

=  3.0 

Input,  as  desired,  the  specific  heat  (Btu/lb-°F) , 

Density  (lb/ft  )  and  emittance 

BSPHT 

=  0.2 

Specific  heat  of  body  skin  material  (Btu/lb-°F), 

default  value  (At) 

BDENS 

=  173.0 

Density  of  body  skin  material  (lb/ft3),  default  value  (At) 

BEPS 

=  0.2 

Infrared  emittance  of  body  skin  material,  default  value 

m 

BTHICK 

=  0.1 

Thickness  of  body  skin  (inch),  default  value 

WMATL 

Wing  skin  material  flag 

=  1.0 

Aluminum,  default  value 
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=  2.0 


Titanium 


=  3.0  Input,  as  desired,  the  specific  heat  (Btu/lb-°F), 

density  (Ib/ft  )  and  emittance 

WSPHT  =  0.2  Specific  heat  of  wing  skin  material  (Btu/lb-°f),  default 

value  (At) 

WDENS  =  173.0  Density  of  wing  skin  material  (lb/ft^),  default  value 

(At) 

WEPS  =  0.2  Infrared  emittance  of  wing  skin  material,  default 

value  (At) 

WTHICK  =  0.1  Thickness  of  wing  skin  (inch),  default  value 

4.2.2  Blunted  Nose 

IN0SE  =  0  Sharp  nose  case 

=  1  Slightly  blunted  nose 

RNOSE  Nose  radius 
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Section  5 
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TEST  CASES 


Test  cases  were  selected  for  aerodynamic  heating  analysis  and  nose  pressure 
modification  using  the  present  code  for  the  purpose  of  verifying  the  theo¬ 
retical  approach  and  the  code  operation. 

5.1  SHARP  CONE  AT  fi  =  4.98 

00 

Results  of  computations  for  heat  transfer  to  a  sharp,  five  degree  semi-apex 
angle  cone  at  =  4.98  are  shown  in  Figure  5-1  compared  with  experimental 
data  from  Reference  14.  The  calculated  results  demonstrate  a  transition 
criterion  of  a  local  Reynolds  number  between  4  x  10^  and  6  x  10^  to 
correspond  with  transition  observed  in  the  test  data.  Comparison  of  calcu¬ 
lated  results  with  experimental  data  for  both  laminar  and  turbulent  boundary 
layer  heat  transfer  is  seen  to  be  very  good. 

5.2  PRESSURE  DISTRIBUTION  OF  A  SLIGHTLY  BLUNTED  NOSE  MISSILE 

For  analyzing  a  missile  with  a  slightly  blunted  nose,  the  test  results  of  a 
circular  cone  with  a  15°  half-cone  angle  at  M  =  1.9  as  reported  in 
Reference  15  was  available.  In  Reference  16,  a  test  for  an  ogive  with  a 
blunted  nose  was  conducted  at  M  =  2.17  and  4.84.  In  Figures  5-2,  5-3  and 
5-4  their  measured  pressure  data  were  compared  with  the  calculated  results 
generated  by  the  code.  It  can  be  seen  that  the  correlation  is  good. 

For  subsonic  speeds,  test  results  for  a  hemisphere-cylinder  at  =  0.8  were 
reported  in  Reference  18.  In  Figure  5-5,  it  can  be  seen  that  calculated  results 
compare  favorably  with  these  data. 


The 

two 


input  requirement 
additional  inputs 


for  the  blunted  nose  code  is  relatively  simple, 
are  needed,  i.e.,  INOSE  and  RNOSE . 
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Figure  5-3.  Pressure  Coefficient  of  a  Slightly  Blunted  Ogive 
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Figure  5-4.  Pressure  Coefficient  of  a  Slightly  Blunted 


'O  zo  So  -<t-o  so 
ANGULAR  DISTANCE  FROM  STAGNATION  POINT  -  DEGREE 


Figure  5-5.  Pressure  Coefficient  of  Hemisphere-Cylinder 


Section  6 

SAMPLE  INPUTS  AND  OUTPUTS 


For  the  convenience  of  the  user,  sample  inputs  end  ootputs  from 
are  presented. 


6.1  SAMPLE  INPUTS 
6  1  1  Aerodynamic  Heating  Case 

For' clarity,  prid  points  end  eights  selected  for  the  sample  test  case  are 
shown  in  the  sketch  (only  a  few  are  shown  for  demonstration  purposes). 


expanded  view  of  nose  elements 


A. 
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6.1.2  Blunted  Nose  Case 

In  the  sketch  shown,  the  blunted  nose  ogive  has  its  elements  and  grid  points 
as  given  by  the  input.  With  the  method  of  analysis  for  the  blunted  nose 
restricted  to  a  slightly  blunted  body,  the  number  of  elements  can  be  small  as 
used  for  the  test  case. 


EXPANDED  VIEW  OF  NOSE  ELEMENTS 
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6.2  OUTPUT  OF  CALCULATIONS 

Printouts  of  the  results  as  calculated  using  the  code  are  reproduced  in  the 
following  pages.  For  the  heating  analysis  case  only  a  few  selected  elements 
were  used  for  demonstration  purposes.  In  application  when  a  large  number  of 
elements  is  involved,  the  temperature  history  will  be  written  on  tape  12  to 
reduce  the  punched  output  and  facilitate  the  handling  of  output. 


6.2.1  Aerodynamic  Heating  Analysis 
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PRESSURES,  FORCES  AND  MOMENTS  ON  ISOLATE  BOOT 

CO  =  .0121  CL  =  .0000  CB  =  0.0000  CPM  =  .0000 


TEMPERATURE  <R)  ON  THE  CONTROL  PANEL 
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6.2.2  Blunted  Nose  Pressure  Distribution 


HMflb 


OQ  «JJK>  •  • 

♦  T  ©-0-00 

yiui  K><\» 

«- »—  ^ 

•  •  •  OO 

%  «  • 

oo°° 

99 

UJIU  IUN  •  • 

r-r-  ift<#00 

•  •  no 

•rvj  *  * 

%  %  »oo 

99 

wut  *  o  *  * 

nr-  Ui  ♦  OO 

«  •  «-ui  •  • 

«(MOO 
%  %  Mn 

r-n  *IO  %  % 

99  Si9°. 

""  8  -00 

•  •  ♦  r—  ^  % 

».  RJ9C?S 

nr*  >uiOO 

.  99 

*  .  <=  .  o'  S£c?° 

OOOO  •■(MO  *00 

It  II  II  N  II  It  II 


O  OO 
O  * 

MO  *  * 
•OOO 


M.  .OO 

n  % 

OO  »  * 
•-OOO 


9899 


r-OOO 
0*0  •  • 

I  »oo 

ui 

ofSS^ 

•*oo 

UI 

VM  %  U 

o§oo 


O  OOOOOC'CCOOOCCC 


O  OO^WXOOCOCl 

o  oooooc^'i'cooo! 


5  *  * 

O  10  ►< 


*  •  %  %  O  '  ^  o  ^  ^  *  %*%%%%*%***%*J 

OOO  ©  ^  o  ©  %  ^  o  o  oooccccccoooocl 

OOO  0000000*^0  COOOO^<COOOQO&| 


ai 


oooooooooooooooooooooooooooo 


oooooooooooooooooooooooooooo 

OOOOOOOOOOOOOOOOOOOOOOOOOOOO 


cr-ooooooooooooo 

0000000*0000000 

OOOOOOCOC'OOOOO 

000000*0000 oooo 


oooooooooooooooooooooooooooo 

oooooooooooooooooooooooooooo 


00000000000000*0  c 


000000*0-0000000 


oooooooooooooo 


oooooooooooooooooooooooooooo 

oooooooooooooooooooooooooooo 

oooooooooooooooooooooooooooo 

oooooooooooooooooooooooooooo 


000000000000000-0000 


OOO'  0000000*0000 


oooooooooooooooooooooooooooo 


oooooooooooooo 


oooooooooooooooooooooooooooo 

oooooooooooooooooooooooooooo 


oooooooooooooooooooooooooooo 


ooooooooooooooooooo- 


0000000000-00000  ooooooooooooooooooooooooooooo  ooooooooooooooooooo 


0000-00000000000  ooooooooooooooooooooooooooooo  OOOOOOOOOC^OOOOOOOO 
OOOOOOOOOOOOOOO  OOOOOOOOOOOOOOOOOOOOOOOOOOOOO  OOOOOOOOOOOOOOOOOOO 


cxrooooooooooooo  ooooooooooooooooooooooooooooo  ooooooooooooooooooo 


OOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOO 


OOOOOOOOOOOOOOO 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

C-’OOOOOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


ooooooooooooooooooo 


ooooooooooooooooooo 

ooooooooooooooooooo- 


ocr-oooooc-ooooooooooooc-oooooooo 

ooooooooooooooooooooooooooooo 


OOOOOC'CCOCO'-OOO 

ooooooooooooooooooooooooooooo 

OOOOOOOOOCOGOC'OOOOCj 

000000*000000000 

ooooooooooooooooooooooooooooo 


00000000000000-0*0000-1 

OOOOOOOOOOOOOOO  00c-00000000000000000000000000  OOOOOOO-  IOOOCKjOCOQOC 


ooooooc^oooooooo  ooooooooooooooooooooooooooooo 
OOOOOOOOOOOOOOO  ooooooooooooooooooooooooooooo 


ooooooooooooooooooo 


63 


WHATL  -  0.0 


BOOT  YAH  *  0.00  0E6 

boot  mu  *  ojjo  m 


PRESSURES,  FORCES  AND  MOMENTS  ON  BODY 


o 


66 


s 


rvjKwDK  fs.  0>  «r 

WQ'J  0^>  «— 
'•O^*^0>'»0*0'»«*nw0«>^» 
»00^0'*0'»,OM<0'M'0»-0*--N 
Ww*-  *-00©0*->  *-  W  W»n*"«  WO 

!••••••«  »**- 


fSWO«^  *-C  OfOO’  r-f-o*  OO  -<*  «#■ 
'»O^*^»0**#C'«0a’0KW5O«»«-« 

•  I  I  «  1  »  •  1  I  r- 


Ww*-  •>OCCO»>  *“WW»0^»WO 

•  !••••••»«- 


0'<><>»OOtA^fO«»K  ■00'*0*T»r- 
O'#  ©*■#  «0*0®>  •Nl'O  «-  ©•—  w 

MMr-  ^OOCO^«-M(MnrMMO 

iia  •  i  r  i*  •*  •  •- 


^OaK^Ult-Kt^PwWO^Q^r 

**"»©'#  O^1  O'*  WtOWM-O^-Or-fJ 
fMI\)»*r*OOOO^^M«V»0K|(\»0 

•  ••••••••*■“ 


f\iK,«)fvr-oOK^O  •— *— O*  ©O  «*-^r 
O'^0kKW>iri*>Ki^K«0rr>*Q>Af- 

*N«NJrf-OOOOf-^<Vir\inrir«0 

I  II  II  •  I  I  l«* 


^^^rt^»T*r-K»«#N®0*rO»rtf 

^r^^OOOOr^NtSn^rgO 

•  I  •  I  (  •  I  I  »’r- 


•VNr-'O/NKp^^irsr^**  N«»»Cv 


Section  7 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  MDAC/NSWC  aerodynamic  code  has  been  extended  to  treat  (a)  aerodynamic  heating 
at  high  supersonic  speeds  and  (b)  nose  pressure  of  a  missile  with  a  slightly 
blunted  nose.  The  input  modifications  have  been  kept  to  a  minimum.  In  the 
aerodynamic  heating  case,  only  a  few  of  the  appropriate  material  properties 
are  required  with  the  default  values  set  for  aluminum.  Only  two  new  inputs 
are  necessary  for  the  blunt  nose  case. 

As  developed,  the  code  represents  a  unique  and  versatile  tool  with  missile 
analysis  capability  that  is  paramount  to  design  studies.  The  advantage  of  the 
output  generation  options  available,  in  the  form  of  NASTRAN  input  specifications 
for  stress  analysis,  needs  no  explanation.  Additionally,  the  small  computing 
effort  required  for  each  case,  must  encourage  its  use. 

When  analyzing  the  operational  capability  of  a  missile,  often  times  only  a  few 
selected  points  on  its  planned  trajectory  are  used  for  study  and  evaluation. 

With  the  advent  of  the  cruise  missile,  considerable  variation  in  the  flight 
environment  may  be  encountered.  To  provide  assessment  of  its  performance  under 
this  situation,  many  cases  must  be  analyzed.  For  this  purpose,  a  flexible 
restart  capability  that  can  meet  all  the  changes  required  for  every  case  while 
eliminating  unnecessary  computational  duplication  will  be  most  welcome.  The 
present  code  is  well  suited  for  modification  to  this  end. 
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